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摘 要 : 以 提取 1990 一 2020 年 植被 覆盖 度 .遥感 4 


E 态 指数 、 人 类 干扰 指数 等 还 感 生态 指标 来 反映 塔 
里 木 河 干流 生态 输 水 工程 实施 前 后 ,生态 状况 变化 趋势 。 将 还 感 生态 指标 作为 驱动 因子 ,提出 一 
种 改进 的 生态 系统 服务 价值 计算 法 ,量化 供给 、 调 节 、 支 持 、 文 化 功能 的 生态 系统 服务 价值 ,分 析 


各 功能 间 权 衡 和 协同 效应 的 动态 演变 关系 ,对 1990 一 2020 年 累积 生态 效益 进行 估算 。 结 果 表 明 : 
(1) 生态 输 水 以 来 ,干流 有 近 1/3 的 区 域 植 被 覆盖 度 增加 , 较 低 、 中 、 较 高 、 高 植被 覆盖 度 面积 占 比 分 
别 增加 17%、5%、2%、2.9%。(2) 上 游 生 态 系 统 服务 价值 先 增加 后 渐 趋 稳定 ,中 下 游 生态 系统 服务 
价值 先 增加 后 减 小 再 增加 ,表明 中 下 游 对 生态 输 水 的 响应 存在 一 段 滞后 期 , 且 受 来 水 水 量 影响 较 
大 。(3) 调节 和 和 支持 功能 之 间 存 在 协同 效应 ,供给 与 调节 和 支持 功能 之 间 存 在 权衡 效应 ,(4) 生态 
输 水 后 期 ,干流 上 、 中 、 下 游标 积 生态 效益 均 表 现 出 “边际 效益 谤 减 规律 ”, 据 此 从 生态 修复 角度 出 
发 ,给 出 断面 年 径流 量 的 适宜 范围 ,上 游 不 超过 42.5x10* m°, 中游 不 超过 21.5x10: m, 下游 不 超过 


3.5x10' m。 研 究 结果 可 为 优化 生态 水 水 量 配置 提供 科学 指导 。 
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工业 化 城市 化 进程 的 迅猛 发 展 导致 人 地 矛盾 
愈演愈烈 ,引发 资源 紧缺 .环境 污染 RE RE 
地 退化 、. 生 物 多 样 性 丧失 等 一 系列 生态 环境 问题 。 
国际 上 , 早 在 20 世 纪 中 叶 已 有 国家 实施 生态 工程 来 
应 对 生态 问题 "1 。20 世 纪 70 年 代 , 我 国 开始 实施 
如 "三 北 "防护 林 工 程 、. 退 耕 还 林 还 草 工程 .生态 输 
水 工程 等 。 目 前 ,生态 修复 已 经 上 升 为 国家 战略 ， 
2017 年 ,党 的 十 九 大 提出 生态 保护 修复 职责 。2020 
年 ,国家 发 展 改 革 委 .自然 资源 部 联合 印发 4 全国 重 
要 生态 系统 保护 和 修复 重大 工程 总 体 规划 (2021 一 
2035 年 )》, 是 从 国家 层面 统筹 生态 保护 与 修复 的 综 
合 性 规划 中 。 截 至 2021 年 全 球 生 态 恢 复 和 保护 面 
积 已 接近 陆地 总 面积 的 17%5。 生 态 修 复 和 保护 工 
程 在 一 定 程 度 上 实现 了 污染 治理 土壤 保持 、 生 物 
多 样 性 恢复 等 生态 目标 。 但 有 学 者 对 全 球 两 百 多 个 
生态 恢复 案例 分 析 发 现 , 仅 有 1/3 的 生态 系统 得 到 了 
有 效 恢 复 , 仍 有 许多 生态 工程 未 达到 预期 效果 “。 
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对 生态 工程 进行 效益 评估 ,能 够 全 面 及 时 地 把 握 生 
态 修 复 成 效 及 存在 的 问题 ,有 助 于 适时 为 管理 决策 
者 提供 调整 策略 。 

生态 工程 能 否 取 得 预期 成 效 , 很 大 程度 上 取决 
于 资源 投入 与 开发 强度 的 适宜 性 ””。 资 源 投入 产 
生效 益 的 多 少 往 往 与 生态 工程 所 处 的 开发 阶段 有 
关 , 根 据 经 济 学 中 “边际 效益 递减 规律 "” ,生态 工程 
从 实施 初期 至 后 期 ,每 增加 单位 立方 米 水 资源 所 产 
生 的 生态 增 量 效 益 呈 递减 趋势 。 雇 淑敏 等 "对 塔 
里 木 河 流域 生态 输 水 的 累积 生态 响应 进行 研究 , 结 
果 表 明 地 下 水 埋 深 和 植被 恢复 随 累 积 生态 输 水 量 
的 增长 性 逐渐 减弱 , 且 具 有 一 定 滞后 性 。 何 理 等 " 
对 “ 引 江 济 人 梨 ” 水 系 连 通 工程 综合 效应 进行 评价 , 结 
果 表 明 连 通 效 果 遵 循 “ 边 际 效益 递减 规律 " 且 具 沛 
后 性 。Ling 等 “运用 这 一 规律 识别 了 干旱 缺 水 地 区 
水 土 资源 开发 的 临界 值 。 对 生态 工程 的 累积 效益 
进行 计算 有 助 于 识别 生态 工程 的 实施 阶段 ,适时 调 
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整 资源 投入 与 开发 策略 ,以 实 
寺 续 发 展 。 

广义 的 生态 效益 指 的 是 生态 系统 对 人 类 生产 
生活 和 环境 产生 的 有 益 影响 和 有 利 效果 外 ,其 内 
涵 是 指 生 态 系统 中 水 、 土 、 气 、 生 等 生态 要 素 通 过 能 
量 流动 .物质 循环 和 信息 传递 等 生态 过 程 ,提供 各 
类 生态 系统 服务 功能 对 人 类 社会 和 自然 环境 产生 
的 积极 作用 。 本 研究 中 生态 效益 通过 生态 系统 服 
务 价值 进行 定量 计算 。 生 态 系 统 服务 价值 计算 法 
主要 有 基于 市 场 理论 的 价值 评估 法 ,基于 生态 模型 
和 当量 因子 的 物质 转换 法 ,能 值 转换 法 等 学 。 目 前 
塔里木 河 生态 系统 服务 价值 的 相关 研究 大 多 运用 
当量 因子 法 , 张 帅 等 “探讨 了 生态 输 水 对 土地 利用 
以 及 生态 系统 服务 价值 的 影响 。 音 静 等 "揭示 了 
土地 利用 变化 对 生态 系统 服务 价值 的 内 在 影响 机 
制 ,提出 了 各 类 土地 利用 最 优 空间 占 比 格局 。 白 元 
等 引 分 析 了 土地 利用 结构 变化 对 生态 系统 服务 功 
能 及 价值 结构 的 影响 ,识别 了 生态 系统 服务 价值 对 
不 同 土地 利用 的 敏感 性 。 然 而 当量 因子 法 在 计算 
生态 系统 服务 价值 时 ,忽略 了 同 种 土地 利用 内 生态 
系统 服务 价值 的 空间 异 质 性 ,因而 有 学 者 对 其 进行 
了 改进 。Ling 等 所 结合 干旱 区 生态 系统 的 特征 , 选 
取 植 被 覆盖 度 、 温 度 植被 干旱 指数 、 人 类 干扰 指数 
作为 生态 系统 服务 价值 计算 的 驱动 因子 ,增强 了 生 
态 系统 服务 价值 的 空间 异 质 性 表达 。 在 此 基础 上 ， 
本 文 进一步 将 其 改进 , 除 考 虑 将 植被 覆盖 和 人 类 干 
扰 指数 作为 驱动 因子 外 ,引入 了 遥感 生态 指数 ,是 
从 热度 . 绿 度 FE .湿度 4 个 方面 对 生态 状况 的 综 
合 反映 ,充分 表征 水 土气、 热 `. 生 等 生态 因子 对 生 
态 系统 服务 价值 的 影响 ,提高 生态 系统 服务 价值 对 
土地 利用 变化 响应 的 敏感 性 。 

20 世 纪 70 年 代 , 塔 里 木 河 下 游 断 流 , 河 岸 植 被 
退化 ,面临 土地 荒漠 化 威胁 。2000 年 起 ,政府 部 门 
开展 了 以 生态 恢复 为 目的 的 生态 输 水 工程 下 。 截 
止 到 2022 年 ,已 实施 生态 输 水 23 次 ,塔里木 河 干流 
生态 发 生 了 极 大 改变 。 本 研究 运用 长 时 间 序 列 
(1990 一 2020 年 ) 遥 感 生态 指标 ,反映 塔里木 河 干流 
生态 输 水 工程 实施 前 后 生态 变化 状况 ; 运用 改进 的 
生态 系统 服务 价值 计算 法 衡量 了 生态 输 水 工程 实 
施 前 后 ,供给 、 调 节 、 支 持 、 文 化 四 类 生态 系统 服务 
功能 之 间 的 权衡 与 协同 关系 的 动态 变化 ;通过 计算 
30 a 来 累积 生态 效益 的 变化 趋势 ,为 今后 生态 输 水 


现 资源 高 效 利用 与 可 


决策 提供 依据 。 
1 数据 与 方法 


11 研究 区 概况 
研究 区 域 为 塔里木 河 干流 (39?30' ~43°08' N, 
73°10'~94°05'E) ,全 长 1321 km™” , 英 巴 扎 、 恰 拉 分 
别 为 中 下游 分 界 点 (图 1)。 塔 里 木 河 深 居 内 陆 , 降 
水 稀少 ,蒸发 旺盛 ,是 典型 的 温带 大 陆 性 干旱 气 
候 。 水 资源 短缺 ,生态 环境 极为 脆弱 。 河 流 两 岸 分 
布 有 广阔 的 荒漠 河岸 林 , 植 被 类 型 主要 有 胡杨 (Pop- 
ulus euphratica) , 1% #9) (Tamarix chinensis) 、 EA 4 hil 
(Halimodendron halodendron) , P33 (Phragmites aus- 
tralis ) 37 Ta PR (Apocynum venetum) IKE H FE (Glyc- 
yrrhiza inflata) FP. 1970s H ERW F WF TA, I 
HARER GEES Tb ie, V T E EER , 
生态 一 旦 退化 ,下 游 将 由 绿洲 变 为 沙漠 ,出 现 不 可 
逆转 的 局 面 。 生 态 输 水 工程 自 2000 年 开始 实施 , 主 
要 指 通过 流域 水 资源 统一 管理 ,增加 各 源流 汇 入 干 
流 的 水 量 , 在 多 年 平均 来 水 条 件 下 ,塔里木 河 干流 
阿拉 尔 断 面 来 水 量 达到 46.5x10 mi, 大 西海 子 断 面 
下 泄 水量 3.5x10* m ,水流 到 台 特 玛 湖 。 生 态 输 水 工 
程 的 目的 是 保护 和 恢复 塔里木 河 干流 上 中 游 林 草 植 
被 ,恢复 下 游 濒危 生态 汪 。 多 年 以 来 生态 输 水 缺乏 
精细 化 调控 方案 , 输 水 方式 随机 日 粗放 ”。 
1.2 数据 来 源 与 处 理 
研究 数据 主要 来 源 于 帝 感 数据 和 实测 数据 
(1990 一 2020 年 )。 实 测 数据 包括 阿拉 尔 新 其 满 、 
英 巴 扎 、 乌 斯 满 、 恰 拉 水 文 站 月 径流 数据 。 融 感 数 
据 被 广泛 应 用 于 评估 大 尺度 、 长 时 间 序 列 、 周 期 性 
的 生态 演变 。 遥 感 数据 包括 植被 覆盖 度 (Fraction 
of vegetation coverage,FVC ) ， 遥感 生态 指数 (Remote 
sensing ecological index, RSEI) 和 人 类 干扰 指数 (Hu- 
man disturbance index, HDI). FVC 和 RSEI 均 来 源 
于 Landsat 影 像 ,空间 分 状 率 为 30 m。HDI 提 取 自 逐 
年 土地 禾 被 产品 [Landsat-derived annual land cover 
product of China (CLCD), https://zenodo.org/record/ 
5816591]” ,空间 分 辩 率 为 30 m。 
FVC 可 以 反应 植被 长 势 和 空间 分 布 ,通过 归 一 
化 植被 指数 (Normalized difference vegetation index, 
NDVIT) 计 算得 出 : 
FVC=(NDVI— NDVI,)/NDVI,..- NDVI,,) (1) 
NDVI=(NIR ~ Red)/(NIR + Red) (2) 
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(b) 阿拉 尔 年 径流 量 及 时 段 平均 值 
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(d) 恰 拉 年 径流 量 及 时 段 平 均值 
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图 1 研究 区 概况 及 断面 年 径流 量 


Fig.1 Sketch map of the study area and annual runoff of the sections 


式 中 :NDVI 为 地 表 纯 裸 地 (或 无 植被 覆盖 ) 时 的 最 
小 NDVI 值 ;NDVI 为 地 表 完 全 植被 覆盖 时 的 最 大 
NDVI {Ë ; NIR X Landsat 影像 近 红 外 波段 (0.76-~ 
0.96 hmy) 的 反射 率 ;Red 为 Landsat 影 像 红 色 波 段 的 反 
射 率 , 波 长 0.62~0.69 hm。NDVI 大 于 0 分 类 为 植被 。 
FVC 可 分 为 5 个 等 级 : 低 [0.0, 0.2) , 较 低 [0.2, 0.4) ,中 等 
[0.4, 0.6) , 较 高 [0.6, 0.8) ,高 [0.8, 1.0]。 
RSET 是 热度 . 绿 度 FE .湿度 指数 经 主 成 分 分 
析 后 得 到 的 结果 ,可 以 从 多 个 方面 反映 生态 状况 ， 
避免 了 指标 之 间 的 相关 性 以 及 人 为 赋 权 的 主观 性 ， 
多 被 用 于 生态 质量 评估 。 通 过 将 多 个 变量 的 特征 
光谱 空间 坐标 轴 进 行 旋转 ,去 掉 各 指标 之 间 的 相关 
性 ,把 主要 信息 集中 于 1~2 个 主 成 分 上 ,可 以 实现 少 
数 重要 变量 的 多 维 数据 压缩 于 。 
RSEI=/(,G,D, W) (3) 
RSEI = f(LST, VI, NDBSI, Wet) (4) 
式 中 :7 为 热度 ;C 为 绿 度 ;也 为 干 度 ; 允 为 湿度 ;LST 
为 地 表 温 度 ; VI 为 植被 指数 ( 同 NDVI) ; NDBSI 为 干 
度 指数 ; Wet 为 湿度 分 量 。 
LST = T/[1 +(ATip)Ine] (5) 
T= K,/\n(K,/L, + 1) (6) 


L,=gain X DN + bias (7) 
式 中 :LST 为 热度 指数 ;4 为 Landsat-7 ETM+ 影 像 6 波 
段 的 中 心 波长 (4=11.45 um) pN 1.438107 m- 
K;e 为 地 表 比 辐射 率 , 不 同 地 物 取 值 不 同 , 水 体 为 
0.99, 植 被 为 0.96, 其余 取 值 0.92。L 为 Landsat-7 
ETM+ 影 像 6 波 段 的 像 元 在 传感器 处 的 辐射 值 ; DN 
为 像 元 灰 度 值 ;gain 和 bias 为 别 为 6 波段 的 增益 值 与 
偏 置 值 , 可 以 从 影像 头 文件 获得 ;7 为 传感器 处 温度 
1E; Ki All KIA ERAR; K,=606.09 We m+ sr， 
um” , K,=1282.71 K- 

NDBSI =(IBI+ SD/2 (8) 


IBI = {2p,/(p; +p.) -[p./(0. +03) +P2/(p2+p5)]}/ 


(9) 
{2ps/(p; +p,) + [pi/(ps +p;) +ps/(p +p} ? 


SI= [(p; +p;) 一 (ps +p,)|/[(P; +p) + (ps +p,)| (10) 


式 中 : 干 度 指数 (NDBST) 由 建筑 指数 (IBI) 和 土壤 指 
数 (SI) 合 成 ; p; G=1, 2, 3, 4,5, 7) 分 别 为 Landsat-7 
ETM+ 影 像 对 应 波段 的 反射 率 。 

Wet = 0.2626p +0.2141p, + 0.0926p, + 


0.0656p, - 0.7629p, - 0.5388p, (11) 
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式 中 :Wet 为 湿度 指数 。 

RSEI 的 值 在 0~1 之 间 , 越 接近 于 1 表明 生态 状 
况 越 好 ,反之 越 差 。RSEI 可 将 生态 状况 分 为 5 个 等 
级 : 差 [0.0, 0.2) ; 低 [0.2, 0.4) , 中 [0.4, 0.6) , [0.6, 
0.8) , 优 [0.8, 1.0]。 

HDI 用 于 描述 人 类 活动 对 区 域内 不 同 土地 利用 
成 分 的 干扰 强度 ,可 根据 土地 利用 数据 计算 。CLCD 
数据 集 将 土地 利用 分 为 耕地 、 林 地 (乔木 ) 灌木 林 
地 草地 水域. 雪 / 冰 、 裸 地 ,建设 用 地 和 湿地 9 类 。 
HDI 的 计算 公式 为 


HDI= > 4,P,/4 
k=1 


式 中 :m 为 土地 利用 类 型 的 数量 ; A, 为 土地 利用 类 
型 的 总 面积 (km ); P, 为 土地 利用 类 型 的 人 为 干 
扰 强度 参数 (耕地 、 林 地 .草地 .水域 .湿地 和 裸 地 的 
参数 分 别 为 0.54 .0.11.0.24 .0.12 .0.15 和 0.08) ;4 为 
各 类 土地 利用 类 型 的 总 面积 。 结 合 各 类 土地 利用 
的 人 类 干扰 强度 参数 ,将 CLCD 土地 利用 数据 集 重 
分 为 6 类 :耕地 林地 (林地 、 治 木林 )、 草 地 水 域 、 
湿地 和 裸 地 ( 裸 地 、 建 设 用 地 、 雪 / 冰 )。HDI 被 分 为 
4 个 等 级 : 低 [0.0, 0.1) , 较 低 [0.1 0.2) ,中 等 [0.2, 
0.3) ,高 =0.3 ,等 级 越 高 代表 人 类 干扰 强度 越 大 1。 


(12) 


表 1 塔里木 河 干流 各 驱动 因子 生态 系统 服务 功 


1.3 研究 方法 
1.3.1 改进 的 生态 系统 服务 价值 计算 法 改进 的 生 
态 系统 服务 价值 计算 法 ,在 当量 因子 法 计算 的 基础 
EU ,综合 考虑 了 FVC、RSEI 和 HDI 因 子 ,避免 了 
同 种 土地 利用 内 ESV 相同 的 局 限 性 ,有 效 表 达 了 
ESV 的 空间 异 质 信 息 ,增强 了 对 环境 变化 响应 的 敏 
感性 。 本 研究 在 Costanza 等 ”和 白 元 等 中 提出 的 
当量 因子 法 的 基础 上 ,结合 塔里木 河流 域 实际 , 确 
定 了 6 种 土地 利用 .9 种 生态 系统 服务 价值 的 当量 
HFR, 
ESV 的 计算 公式 如 下 : 
ESV(x)= S ESV f) 
式 中 : ESV(x) 为 像 元 x 的 生态 系统 服务 价值 (元 ); 
ESV (x) 为 像 元 x 第 i 种 生态 系统 服务 功能 的 价值 。9 
种 生态 系统 服务 功能 分 别 是 大 气 调 节 、 气 候 调 节 、 
水 源 涵养 .废物 处 理 、 水 土 保持 .生物 多 样 性 维持 、 
食品 供给 、 原 材料 供给 娱乐 和 文化 。 
ESV (x)= YD, (2) x EV,(x) x4, XA 


(13) 


(14) 


式 中 : D(x) 为 第 1 个 驱动 因子 的 第 ;种 生态 系统 服 
务 功 能 的 影响 系数 , 取 值 范围 为 0~1, 各 驱动 因子 的 
影响 系数 见 表 1 中 ,3 种 驱动 因子 分 别 是 FVC 、HDI、 


能 影响 系数 和 权重 


Tab.1 Influence coefficients and weights of ecosystem services for each driving factor in the mainstream of Tarim River 


生态 系统 服务 功能 

驱动 权重 分 级 调节 支持 供给 文化 
因子 KA AR KE 废物 xE 生物 多 样 食物 ”原材料 。 娱乐 和 
调节 调控 涵养 处 理 保持 性 维持 供给 供给 文化 

植被 覆 0364 0.0-0.1 0.01 0.01 0.01 0.05 0.05 0.05 0.05 0.05 0.10 
an JS 01-03 035 035 0.25 0.40 0.40 0.20 0.25 0.40 0.30 
03-05 060 0.60 0.55 0.60 0.60 0.50 0.50 0.60 0.55 

05-07 085 085 0.85 0.80 0.80 0.85 0.75 0.80 0.80 

07-10 1.00 1.00 100 1.00 1.00 1.00 1.00 1.00 1.00 

REF 0.332 0.0-0.1 0.01 001 0.01 0.01 0.05 0.01 0.01 0.01 0.05 
扰 指 数 0.1~0.2 1.00 1.00 1.00 1.00 1.00 1.00 0.50 1.00 1.00 
0.2-0.3 085 085 0.90 0.75 0.80 0.85 0.75 0.55 0.50 

03-04 065 0.65 0.55 0.80 0.75 0.65 0.85 0.60 0.25 

04-10 040 040 0.25 0.75 0.70 0.45 1.00 0.45 0.10 

ÆRE 0304 0.0-0.2 010 0.10 0.05 0.10 0.10 0.10 0.05 0.05 0.10 
态 指数 02-04 035 035 0.20 0.20 0.30 0.30 0.25 0.25 0.25 
0.406 060 0.60 0.50 0.40 0.60 0.60 0.55 0.55 0.50 

0.6-0.8 085 085 0.80 0.60 0.90 0.90 0.90 0.90 0.75 

08-10 1.00 1.00 100 1.00 1.00 1.00 1.00 1.00 1.00 
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RSEI。 EV,(x) 为 各 类 土地 利用 类 型 第 i 种 生态 系统 
服务 功能 的 当量 值 '"; 4 为 第 7 个 驱动 因子 的 权重 
值 ;4 为 李 格 分 辨 率 ,30 mx30 mo 

1.3.2 累积 生态 效益 累积 生态 效益 是 指 在 一 定时 
间 范 围 内 ,由 于 生态 保护 和 恢复 措施 的 实施 ,生态 
系统 各 方面 的 改善 和 效益 逐渐 积累 和 增加 的 总 体 
效果 。 与 生态 效益 相 比 ,累积 生态 效益 关注 的 是 保 
护 和 恢复 措施 对 生态 系统 的 长 期 影响 。 通 过 分 析 
累积 生态 效益 随 累积 断面 年 径流 量 的 变化 趋势 ,可 
以 判断 当前 阶段 资源 投入 与 开发 强度 的 适宜 性 。 
本 文 运 用 生态 系统 服务 价值 与 断面 年 径流 量 之 比 ， 
表征 生态 效益 ,将 生态 效益 逐年 累加 得 到 累积 生态 
效益 ,计算 公式 如 下 : 


„o~ ESV, 
FEB,= 之 WOV (15) 


式 中 : EB! 为 第 ;年 的 累积 生态 效益 ; ESV, 为 第 ;年 
的 生态 系统 服务 价值 ; EWCV ,为 第 ;年 的 断面 年 径 
流量 。 EB' 是 从 第 1a 开 始 累加 到 第 ”年 的 生态 效 
益 之 和 。 


2 结果 与 分 析 


2.1 生态 状况 时 空 变 

本 研究 运用 1990—2020 4# FVC, RSEL, HDI X 
分 析 塔 里 木 河 干流 生态 状况 时 空 变化 趋势 。 将 研 
究 区 域 划分 为 上 游 中游. 下游 3 个 河 段 。 按 照 生 态 
输 水 工程 实施 阶段 ,将 研究 时 段 划分 为 3 个 时 段 : 
1990—1999 4F , 2000—2010 年 和 2011 一 2020 年 分 
别 代表 生态 输 水 工程 实施 之 前 实施 前 期 和 实施 
后 期 。 

塔里木 河和 干 流 FVC 的 结果 显示 (图 2) ,生态 输 
水 工程 实施 之 后 上 .中 .下游 FVC 整 体 提 升 。2011 
年 后 相 较 2000 年 前 干流 上 P. FYHR FVC 面积 
比分 别 减 少 45% 35% 13% ,中 等 以 上 FVC 面积 
比分 别 增加 25%、11%、3%。2000 年 生态 输 水 工程 
实施 以 来 ,干流 有 近 1/3 区 域 FVC 提 升 。 

1990 一 2020 年 RSEI 结 果 显 示 ( 图 3) ,上游 生态 
状况 好 转 ,中 游 表现 出 下 降 趋 势 ,下游 逐步 改善 。 
对 比 1990 一 1999 年 和 2011 一 2020 年 的 RSEI, 上 游 
RSEI“ 良 ” 值 区 增加 了 15%。 中 游 RSEI“ 低 ” 值 区 增 
加 了 10%,“ 中 ” 值 区 减少 了 13%, 生 态 状 况 有 下 降 趋 
势 。 下 游 各 区 间 面 积 占 比 变 化 不 大 ,与 1990—1999 
年 相 比 ,2011 一 2020 年 RSEI“ 良 ” 值 区 增加 2%。 生 


态 输 水 工程 导致 中 游 引水 量 减少 所 ,是 影响 中 游 生 
态 状 况 的 原因 之 一 (图 1)。 

HDI 时空 变 化 趋势 结果 显示 (网 4) , 1990 一 
2020 年 上 游 和 中 游人 类 干扰 在 逐渐 增加 ,下 游人 类 
干扰 基本 稳定 。HDI 数 值 越 大 表明 人 类 活动 干扰 程 
度 越 大 。 与 1990 一 1999 年 相 比 ,2011 一 2020 年 上 游 
和 中 游 HDI 高 值 区 面积 占 比分 别 增加 了 13% 和 5%。 
2.2 生态 系统 服务 价值 计算 

对 干流 上 中、 下 游 1990 一 2020 年 供给 、 调 节 、 
支持 .文化 功能 生态 系统 服务 价值 计算 结果 进行 统 
计 , 结 果 表明 ,调节 和 支持 功能 的 价值 大 于 供给 和 
文化 功能 的 价值 ,调节 和 支持 功能 之 间 存 在 协同 效 
应 ,供给 与 调节 支持 功能 之 间 存 在 权衡 效应 (图 
5)。 随 着 生态 输 水 工程 的 实施 ,水 域 面积 和 植被 盖 
度 增 加 ,有 利于 气候 调节 、 水 源 涵养 .水 土 保持 、 生 
物 多 样 性 等 调节 和 支持 功能 的 发 挥 。 塔 里 木 河 供 
给 功能 主要 依靠 农田 生态 系统 ,农田 面积 的 扩张 会 
导致 林 草 地 和 裸 地 等 的 面积 减少 ,因而 供给 功能 与 
调节 和 支持 功能 之 间 呈 权衡 效 应 。 为 了 尽 可 能 地 
发 挥 各 类 生态 系统 服务 功能 的 效益 ,需要 协调 各 类 
土地 利用 的 面积 。 

总 体 而 言 ,生态 系统 服务 价值 上 游 最 大 ,其 次 
是 中 游 和 下 游 。1990 一 2020 年 上 游 各 类 生态 系统 
服务 价值 呈 先 增加 后 稳定 趋势 。 中 游 1990 一 1999 
年 .2000 一 2010 年 ,2011 一 2020 年 各 类 生态 系统 服 
务 价值 呈 先 增加 ,后 缓慢 减少 ,再 逐渐 递增 趋势 。 
下 游 1990 一 1999 年 生态 系统 服务 价值 变化 不 大 ， 
2000 一 2010 年 呈 先 增加 后 减少 趋势 ,2011 一 2020 年 
逐渐 递增 。 表 明生 态 输 水 工程 实施 之 后 ,上 游 生态 
系统 服务 价值 先 增加 后 渐 趋 稳定 ,中 下 游 生态 系统 
服务 价值 对 生态 输 水 的 响应 存在 一 段 灌 后 期 , 且 受 
来 水 水 量 影 响 较 大 (图 1)。 

2.3 生态 效益 分 析 

图 6a~c 表明 生态 系统 服务 价值 变化 趋势 与 断 
面 年 径流 量变 化 趋势 基本 一 致 ,生态 系统 服务 价值 
随 年 径流 量 同步 增 减 。 上 中游 生态 系统 服务 价值 
先 增加 后 稳定 ,下游 生态 系统 服务 价值 波动 上 升 。 
上 中游 生态 系统 服务 价值 不 会 随 着 生态 输 水 工程 
的 实施 持续 上 升 ,上 升 到 一 定 程 度 之 后 增 速 减 慢 趋 
于 稳定 ,表明 生态 系统 对 生态 输 水 的 响应 存在 靖 
值 。2018 一 2020 年 下 游 生态 输 水 水 量 减少 而 生态 
系统 服务 价值 持续 增加 ,表明 生态 输 水 对 于 生态 恢 
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图 2 塔里木 河 干流 植被 覆盖 度 (FVC) 及 其 面积 占 比 


Fig. 2 Fraction of vegetation coverage and area proportion of the mainstream of Tarim River 
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图 6d~f 散 点 图 表示 的 是 累积 生态 效益 ,柱状 图 
表示 的 是 累积 断面 年 径流 量 。 生 态 输 水 之 后 ,上 、 
中 下 游 累 积 生 态 效 益 递增 速率 (k=2.83 ,k=4.11 ,k= 
8.91) 均 大 于 生态 输 水 之 前 (k=1.66, k=2.58, k= 
8.77) ,表明 生态 输 水 后 流域 生态 恢复 效果 显著 。 生 


态 输 水 后 期 上 .中 下游 累 积 生态 效益 递增 速率 逐渐 
趋 缓 ,而 这 一 时 期 累积 径流 量 递增 速率 加 快 ,表明 随 
着 生态 输 水 的 持续 进行 ,生态 输 水 工程 所 产生 的 生 
态 效益 不 会 持续 递增 ,而 是 存在 一 定 的 阔 值 ”> ， 
这 一 现象 符合 经 济 学 中 的 “边际 效益 递减 规律 ”。 
生态 输 水 进行 到 一 定 阶段 之 后 ,单位 立方 米 水 所 产 
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图 3 塔里木 河 干流 遥感 生态 指数 (RSEI) 及 其 面积 占 比 


Fig. 3 Remote sensing ecological index and area proportion of the mainstream of Tarim River 


生 的 生态 效益 将 不 再 递增 ,因此 可 以 将 输 水 水 量 控 PERELU , MTT SLE AS Sad ic BY ts yE A BK 
制 在 一 定 区 间 内 。 根 据 累积 生态 效益 转折 点 出 现 。 恢复 “ 提 质 增 量 ”的 目的 ?31。 

的 时 间 ,对 断面 年 径流 量 的 适宜 范围 建议 如 下 ,上 

游 不 超过 42.5x10* m ,中游 不 超过 21.5xl0m 下游 3 结论 

不 超过 3.5x10: ms。 与 此 同时 ,在 输 水 总 量 一 定时 ， 

充分 保证 植被 萌芽 期 的 生态 需 水 ,结合 不 同 生 态 功 (1) 本 文 基于 1990—2020 年 遥感 影像 提取 
能 区 的 恢复 目标 ,确定 不 同 来 水 频率 下 的 生态 供水 。 FVC、RSEI、HDI, 分 析 塔 里 木 河 干流 生态 状况 时 空 
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4 塔里木 河 干流 人 类 干扰 指数 (HD 了 DD) 及 面积 占 比 


Fig. 4 Human disturbance index and area proportion of the mainstream of Tarim River 


变化 趋势 ;将 以 上 3 个 指标 作为 驱动 因子 提出 了 一 
种 改进 的 生态 系统 服务 价值 计算 法 ,定量 解析 了 
上 中、 下游 生 态 输 水 前 后 4 种 生态 系统 服务 功能 间 
权衡 和 协同 效应 演变 关系 ;分 别 计算 各 河 段 
累积 生态 效益 ,通过 分 析 变 化 趋势 ,对 生态 输 水 实 
践 提 供 指导 。 


(2) 生态 输 水 以 来 干流 有 近 1/3 的 区 域 FVC 提 
升 ,上 游 FVC 提 升 面积 占 比 最 大 ,依次 是 中 游 和 下 
游 。RSEI 结 果 表 明 干 流 上 游 和 下 游 生态 状况 逐渐 
好 转 , 中 游 有 下 降 趋势 。HDI 结 果 显 示 1990—2020 
年 上 游 和 中 游人 类 干扰 逐渐 增加 ,下 游 基 本 稳定 。 

(3) 干流 生态 系统 服务 功能 以 调节 和 支持 功能 
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图 6 塔里木 河 干流 生态 效益 与 径流 量 
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Fig.6 Ecological benefits and runoff of the mainstream of Tarim River 
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占 主导 ,调节 和 支持 功能 存在 协同 效应 ,供给 与 调 
节 和 支持 功能 之 间 存 在 权衡 效应 。 研 究 发 现 生态 
效益 遵循 “边际 效益 递减 规律 ”, 当 累积 断面 年 径流 
量 达到 某 一 阔 值 后 ,无 法 再 通过 增加 输 水 水 量 的 方 
式 获 取 更 多 的 生态 效益 ,此 时 可 将 上 游 中游. 下游 
断面 年 径流 量 维持 在 一 定 范 围 , 上 游 不 超过 42.5x 
10° ma ,中游 不 超过 21.5x10: m ,下 游 不 超过 3.5x 
10: mw?。 在 今后 的 研究 中 可 以 依据 生态 需 水 关键 期 
和 不 同 生态 功能 区 的 恢复 目标 ,确定 不 同 来 水 频率 
下 的 生态 供水 轮 治 规则 。 

(4) 改进 的 生态 系统 服务 价值 计算 法 包含 了 
FVC、RSEI、HDI 3 个 驱动 因子 , 受 数 据 可 获得 性 的 
限制 ,没有 将 土壤 售 水 量 、 地 下 水 埋 深 等 表征 生态 
系统 变化 过 程 的 因子 考虑 在 内 ,可 结合 研究 区 域 特 
征 与 研究 需求 进一步 将 其 改进 。 
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Ecosystem variation and ecological benefits analysis of the 


mainstream of Tarim River 


HU Ruiyuan, CHANG Jianxia, GUO Aijun, WANG Yimin 
(State Key Laboratory Base of Eco-hydraulic in Arid Area, Xi’ an University of Technology, Xi an 710048, Shaanxi, China) 


Abstract: In this study, remote sensing ecological indices such as vegetation coverage, remote sensing ecological 
index, and human disturbance index were extracted from 1990 to 2020 to reflect the evolution trend of ecological 
status before and after the implementation of the ecological water conveyance project of Tarim River. Taking re- 
mote sensing ecological indices as driving factors, an improved ecosystem service value calculation algorithm 
was proposed to quantify the ecosystem service value of provisioning, regulating, supporting, and cultural func- 
tions and to analyze the dynamic evolution relationships between trade-offs and synergistic effects among the 
functions. Cumulative ecological benefits from 1990 to 2020 were estimated. The results showed that the ecologi- 
cal condition of approximately 1/3 of the mainstream of the Tarim River had improved because of the implemen- 
tation of the ecological water conveyance project. The proportion of lower, medium, higher, and high vegetation 
coverage increased by 17%, 5%, 2%, and 2.9%, respectively. The results of the ecosystem service value calcula- 
tion showed that the ecosystem service value of the upstream increased and then stabilized, whereas that of the 
midstream and downstream increased, decreased, and then increased again, indicating a period of lag in the mid- 
stream and downstream in response to the ecological water conveyance and was significantly affected by the 
amount of incoming water. There was a synergistic effect between the regulating and supporting functions and a 
trade-off effect between the provisioning function and the regulating and supporting functions. Cumulative eco- 
logical benefits after 2010 followed the “marginal benefit decreasing law”. Accordingly, for ecological restora- 
tion, the annual runoff should not exceed 42.5x 10° m upstream, 21.5x 10° m midstream, and 3.5 10° m° downstream. 
This study can provide scientific guidance for optimizing the allocation of ecological water quantity. 


Key words: ecological water conveyance project; ecosystem service; ecological benefits; Tarim River 


